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Abstract. It is supposed that the exponential multiplier in the method of the 
non-equilibrium statistical operator (Zubarev's approach) can be considered as a 
distribution density of the past lifetime of the system, and can be replaced by an 
arbitrary distribution function. To specify this distribution the method of maximum 
entropy principle as in [Schonfeldt J-H, Jiminez N, Plastino A R, Plastino A, Casas M 
2007 Physica A 374 573] is used. The obtained distribution is close to exponential one. 
Another approach to the maximum entropy principle, as in [Van der Straeten E and 
Beck C 2008 Phys. Rev. E 78 051101], except exponential distributions yields power- 
like, log-normal distributions, as well as distributions of other kind and transitions 
between them. 
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1. Introduction 

Among possible approaches to the description of non-equilibrium systems the Non- 
equilibrium Statistical Operator Method (NESOM) especially demonstrated its efficiency 
[H El |3j . NESOM provides a very promising technique that implies in a far-reaching 
generalization the statistical methods developed by Boltzmann and Gibbs. NESOM was 
initially built on intuitive and heuristic arguments, apparently it can be incorporated 
within an interesting approach to the rationalization of statistical mechanics, as 
contained in the maximization of (informational statistical) entropy (MaxEnt for short) 
and Bayesian methods. NESOM appears as a very powerful, concise, based on sound 
principles, and elegant formalism of a broad scope to deal with systems arbitrarily far 
from equilibrium. The non-equilibrium statistical operator (NSO) introduced in [H [2J, [3] 
has a form 

POO 

\ng{t)= p qzub {u) In g q (t-u, -u)du, In g q (t, 0) = -^{t)-^ F n (t)P n ; 
Ing^h) = e^ H ^\ng q (t,0)e^ ih h = In Spexp{- ^ F n (t)P n } , 

n 

where 

p q (u) =p qzu b{u) = ee~ £U , u = t-t , (2) 

H is Hamiltonian, \np(t) is the logarithm of the NSO, In p q (ti,t2) is the logarithm 
of the quasi-equilibrium (or relevant) distribution; the first time argument indicates 
the time dependence of the values of the thermodynamic parameters F n ; the second 
time argument t<i in lnp g (ti,t 2 ) denotes the time dependence through the Heizenberg 
representation for dynamical variables P n on which In p q (t, 0) can depend [H [2j [3]. 

In [U El [3] p q {u) = p qzu b{u) = e exp{ — eu}; after the thermodynamic limiting 
transition iV — > oo, V — > oo, N/V = const, e — > 0. From the complete group of 
solutions of Liouville equation (symmetric in time) the subset of retarded "unilateral" 
in time solutions is selected by means of introducing a source K in the Liouville equation 

, , , w r . T . .„ . r r TT i \9H dg dH dg 1 

for lnp(t) (L is Liouville operator; iL = —{si, g\ = 2_^ k — — — — — — — ; pk, and 

\_op k oq k oq k op k \ 

qk are pulses and coordinates of particles; {. . .} are Poisson brackets) 

dln £® +iL\np{t) = -e(ln p(t) - lnp,(t,0)) = K zub . (3) 

In (U [S] a convenient redefinition of the source term is proposed. Although 
infinitesimally small, the source term introduced by Zubarev into the Liouville equation 
is shown to influence the macroscopic behaviour of the system in the sense that the 
corresponding evolution equations do not coincide exactly with those obtained from 
an initial-value problem which corresponds to a definite experimental situation and a 
physical set of macroobservables. 

In [6l [7] it was noted that in place of the function p q {u) = p qzu b{u) = ee~ £U in 
([T]) arbitrary (but having certain properties [S]) weigth functions w(t,t ) can be used. 
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Zubarev's nonequilibrium statistical operator does satisfy Liouville equation, but it 
must be borne in mind that the group of its solutions is composed of two subsets, 
one corresponding to the retarded and second one to the advanced solutions. The 
presence of the weight function w(t,t ) (Abel's kernel (jSJ) in Zubarev's approach) in the 
time-smoothing or quasi-average procedure that has been introduced selects the subset 
of retarded solutions from the total group of solutions of the Liouville equation. This 
consideration is related to the question: how to obtain an irreversible behavior in the 
evolution of the macroscopic state of the system? In the MaxEnt-NESOM approach the 
irreversibility is incorporated from the outset using an ad hoc non-mechanical hypothesis. 
MaxEnt-NESOM yields information on the macrostate of the system at time t, when 
a measurement is performed, including the evolutionary history (in the interval from 
the initial time of preparation to up to time t) by which the system came into that 
state (which introduced a generalization of Kirkwood's time-smoothing formalism [9]). 
Functions w(t, t ) are typically kernels [U [2J [3l [8] that appear in the mathematical theory 
of convergence of integrals. In pill Qj], [12] other interpretation of the functions w(t,t ), 
denoted as p q (u), is given. With the change of function w(t,t ) = p q (u) the form of 
source (j3J) in the Liouville equation also changes. For an arbitrary function p q {u) it 
looks like (gj. 

In [101 [HJ [12] it was noted that the function p q (u) = p qzu b( u ) = ee~ eu in 
NESOM [U El [3] for the non-equilibrium distribution function can be interpreted as 
the exponential probability distribution of the lifetime T of a system. T is a random 
variable of lifetime (time span) from the moment to of its birth till the current moment t; 
e^ 1 = (t — 1 ) = (r), where (T) = J up q {u)du is the average lifetime of the system. This 
time period can be called the time period of getting information about system from its 
past. Instead of the exponential distribution p qzu b{u) ([2]) in ([1]) used in [HO [3] any other 
sample distribution p q (u) could be taken; integration by parts in time is performed at 
/ p q (u)du\ u= o = -1; / pqiujdu^oo = 0. If p q (u) = p qzu b{u) = ee~ £U by (j2J), e = 1/(T) 
the expression for NSO passes in flTJ from [H [2J, [3] . 

The same interpretation of the distribution p q (u) is given in [2], where this value 
is understood as the distribution of the initial moment of time t . Since the random 
(past) lifetime is equal to T = u = t — to, the distribution of the past lifetime u coincides 
with the distribution of the initial time values to- The moment to will be the moment 
of the first passage in the inverse time, if the moment t is taken as initial. In [2] 
the uniform distribution for an initial moment to is chosen, which after the transition 
from Abel integration to Cesaro integration passes to the exponential distribution 
p q (u) = Pqzubiu) = ee~ eu . Such distribution serves as the limiting distribution of the 
lifetime [L3], the first-passage time of a certain level. In the general case it is possible 
to choose a lot of functions for the obvious type of distribution p q (u), which was noted 

in HD1HHH2J. 

In [13] the lifetimes of the system are introduced as random moments of the first- 
passage time till the moment when a random process describing system reaches a certain 
limit, for example, a zero value. In [13J approximate exponential expressions for the 
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probability density function (with a single parameter) and probability distribution of 
lifetime are obtained, and the accuracy of these expressions is estimated. 

In p3] it was noted, that the role of the form of the source term in the Liouville 
equation in the NSO method has never been investigated. In [15] it is stated that 
the exponential distribution is the only one which possesses the Markovian property 
of the absence of afteraction, that is whatever is the actual age of a system, the 
remaining time does not depend on the past and has the same distribution as the 
lifetime itself. It is known [JJ EJ [3] that the Liouville equation for NSO contains the 
source K = K zuh = — e[\np(t) — \np g (t,0)] (J3]) which becomes vanishingly small after 
taking the thermodynamic limit and setting e — > 0, which in the spirit of pU] corresponds 
to the infinitely large lifetime value of an infinitely large system. For a system with finite 
size this source is not equal to zero. In [8] this term enters the modified Liouville operator 
and coincides with the form of Liouville equation suggested by Prigogine [16] (the 
Boltzmann-Prigogine symmetry), when the irreversibility is introduced in the theory 
at the microscopic level. 

In [10] a new interpretation of the method of the NSO is given, in which the 
operation of taking the invariant part [U [2] or the use of an auxiliary "weight function" 
(in the terminology of P [7J [8]) in NSO are treated as averaging the quasi-equilibrium 
statistical operator over the distribution of past lifetime of a system. This approach 
agrees with the approach of the general theory of random processes, the renewal theory, 
and also with the conception of Zubarev work p] where the NSO is conceived as some 
averaging over the initial moment of time. 

The statistical operator depends on the information-gathering interval (t ,t), but 
it must be borne in mind that this is the formal point consisting in that (as Kirkwood 
pointed out) that the description to be built must contain all the previous history 
in the development of the macrostate of the system. In [61 [8] several basic steps 
for the construction of the NESOM formalism are indicated: a third basic step has 
just been introduced, namely, the inclusion of the past history (other terms used 
are retro-effects or historicity) of the macrostate of the dissipative system. A fourth 
basic step needs now to be considered, which is a generalization of Kirkwood's time- 
smoothing procedure: the one that accounts for the past history and future dissipative 
evolution. The time-smoothing procedure introduces a kind of Prigogine's dynamical 
condition for dissipativity. The procedure introduces a kind of evanescent history as 
the system macrostate evolves toward future from the initial condition at time to- The 
function w(t] t ) [HI [7J [8] introduces the time-smoothing procedure. In principle, any 
kernel provided by the mathematical theory of convergence of trigonometrical series 
and transform integrals provides is acceptable for these purposes. Kirkwood, Green, 
Mori [9j [T71 [18] and others have chosen what in mathematical terminology is Fejer (or 
Cesaro-1) kernel. Meanwhile Zubarev introduced the one consisting in Abel's kernel for 
w in Eq. ([lj - which apparently appears to be the best choice, either mathematically but 
mostly physically: that is, taking w(t;to) = ee e ^ 0_i \ where e is a positive infinitesimal 
value which tends to zero after the calculation of averages has been performed, and with 
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to going to — oo. Therefore a process with fading memory is introduced. In Zubarev's 
approach this fading process occurs in an adiabatic-like form towards the remote past: as 
time evolves memory decays exponentially with lifetime e~ l [8]. The approach suggested 
in [10j[TTJ[T2] an d in the present work enables to use a family of functions w(t, t ) = p q {u) 
and makes clear both their physical sense and those physical situations in which one or 
another function w(t,t ) = p q (u) can be used. 

Besides the Zubarev's form of NSO [HOE], the NSO formulation in the Green- 
Mori form [T71 [IE] is known, where one assumes the auxiliary weight function [HI [8] 
to be equal to W(t,t ) = 1 — (t — t )/r; w(t,t ) = dW(t,t )/dto = 1/r; r = t — to- 
After averaging one sets r — > oo. This choice at p q (u = t — to) = w(t, to) coincides with 
the uniform lifetime distribution. The source in the Liouville equation takes the form 
K = In p q /r. In pQ this form of NSO is compared to the Zubarev's form. One could 
name many examples of explicitely setting the function p q (u). Each and every definition 
implies some specific form of the source term K in the Liouville equation, some specific 
form of the modified Liouville operator and NSO [TQl Qj], [12]. Thus the whole family of 
NSO is defined. 

It is possible to make different assumptions about the form of the function of 
p q {u), getting different expressions for the source in the Liouville equation and for non- 
equilibrium characteristics of the system. It is possible to show [12] that certain choices 
of the function of p q (u) result in the changes in non-equilibrium characteristics in the 
limit of infinitely large average lifetimes as well. In [101 [TH [12] an analogy is traced 
to the passage to the thermodynamic limit of systems of infinite size. So explicit form 
of the function p q (u) is important for describing non-equilibrium systems by the NSO 
method. 

Setting the form of the function p q (u) reflects not only the internal properties of 
a system, but also the influence of the environment on an open system, the particular 
character of its interaction with the environment [8]. In [2] a physical interpretation of 
the exponential distribution for the function p q (u) is given: a system evolves freely like an 
isolated system governed by the Liouville operator. Besides that the system undergoes 
random transitions, and the phase point representing the system switches from one 
trajectory to another one with an exponential probability under the influence of the 
"thermostat"; the average intervals between successive push events increase infinitely. 
This takes place if the parameter of the exponential distribution tends to after the 
transition to the thermodynamic limit. Real physical systems have finite sizes. The 
exponential distribution describes completely random systems. The influence of the 
environment on a system can have organized character as well, for example, this is the 
case of systems in a stationary non-equilibrium state with input and output fluxes. The 
character of the interaction with the environment can also vary; therefore different forms 
of the function p q (u) can be used. 

The adequate choice of the function p q (u) is important for correct description of 
the non-equilibrium properties of statistical systems. To find the type of function p q (u), 
it is necessary to resort to some general principles, such as MaxEnt principle. In this 
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work two variants of MaxEnt are used for this purpose: the one introduced in 
for the Liouville equation with a source (Section 2) and that suggested in [26J for the 
superstatistics (Section 3). 

2. Maximum entropy principle for Liouville equations with source 

In this paper we apply the maximum entropy principle for the determination of the 
function p q {u). The same approach was applied in [19] for the evolution equations with 
source terms. In [3, [lOj [TT], [12] a general form for the source in the Liouville equation 
for In p(t) ([3]) is obtained. For our case the source term has the following form 

K = p, (0) In Pq (t, 0) + In p q {t - u, -u)du. (4) 

In [TS] in the Liouville equation the distribution function p(z,t) is written in a form 

N M 

p{-?, t) = Nf ME {!?, t) = — exp{- ^ M*}, (5) 

i=i 

where Ai(z) are M appropriate quantities that are functions of the phase space point 
z\ the quantities A\{z) correspond to the values Pi from The partition function Z 
is given by 

» M 

Z= exp{-J2\iAi}d N z. (6) 

J i=l 

The function / ME (f,i) is normalized to unity: 



f ME (l?,t)d"z = l; (7) 

J p{-t,t)d N z = N(t); ^ = j Kd N z; (8) 

% + UVp = ^ = K- = 0. 

ot dt 

The probability distribution fME(z,t) is the one that maximizes the entropy S[f] 
under the constraints imposed by normalization and relevant mean values (Aj) = 
/ Aipd N z (or a, = (A)i/N). The re-scaled mean values a, and the associated Lagrange 
multipliers Aj are related by the Jayne's relations [201 ETJ 

S= - I f In fd N z = In Z + J2X iai . 

J i 

If we choose for In p(z, t) the function 

poo pt 

\np = \no(t)= I p q (u)\n g q (t—u, —u)du = I p q (t— 1 ) In p q (t ,t —t)dt { 

Jo J -oo 
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from (pQ), which is included in the Liouville equation (without integrating over time) 
with p qzu b — - > p q , K zub — > K, and, following [I], choose 



then 



where 



\i = Pq {t-t )F i {t ), (11) 
lnp(z,t) = p q (t- t ) In p q (t , t -i)dt = (12) 

J — oo 

- [ (J2 X ^ + P«(* - £ o) 111 Z l) rft = 111 f'ME + 111 N > 
J-oo i 

In /j^ = - / Aj Ajrfto - In Z A ; 

J-OO j 

InN = AZ = \nZ x - f p q {t - t ) In Zidt ; (13) 

•/ — oo 

Z A = / exp[- I ^X^dt^af z, Z x = f expf-J^A]^ 

J J-oo i 7 i 

(Fi are taken from (pQ)). The values Z A and in the terminology of [1] are related to the 
partition functions for a non-equilibrium and relevant statistical operator accordingly. 

In [19] for the Liouville equation of the kind (j3J) with constant sources equation one 
gets for dXi/dt 

dX, ,A „ _ 19 



i=i * ^ 

where the Zubarev-Peletminskiy selection rule [221 [2H [H E] 



dt 

i=i 



M 



wVA^J^C^, (i = l,...,M), § = ^ + V(p^)=if (14) 



is used; z,j = 1,2...; the C^- are c-numbers. In other representations the quantities A+ 
can depend on the space variable, that is, when considering local densities of dynamical 
variables, and then the can depend on the space variable as well or be differential 
operators. 

If more complex shape of the source (Hj) is considered, the equation for dXi/dt takes 
on the form 
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We replace the operators — — and — - taking into account fl9l)- f|T3l) by the functional 
differentiation of the kind 

_d 5_ 5 _d 5_ 

which takes off the integration over time. For example 



k 



The relations ([ME]) and (jUJ) are thus hold. If to take into account that J . . . p q d 
J . . . pd N z = (...) in the NSO method, then 



z 



A»=i-W=0. (16) 
dca dai yNj K ' 

Let us consider the integrals in the rhs of ( |T5l) of the form j KB(z)dz, B being an 
arbitrary function of the dynamic variables z, and the source term K taken from (j3J); 

for Eq. USD, dUD K = \p q (0) lnp q {t, 0) + J ^hSLlM \n Pq (t ,t - t)dt ]p. Assume 

that p q {u) does not depend on z. Integrating by parts and assuming p g (u)„^oo - > 0, we 
get: 

/r°° d 
KBd N z = -J p q (u) — {Bln Pq (t-u,-u))du; u = t - 1 . (17) 

d & 

Taking into account f|T6l) and the fact that the operation — > - — eliminates the 

ddi oa,i 

integration by time, the equation (fT5l) takes on the form 

Fi{h) dPq ^~ to) = -p q (t - t )Ci - Pq (t - t )n - p\(t - t )r 2 , (18) 
where Q = X] ; <",,/•}(/,)). 



n = -qIXXJ^p^) ~ lnA ln P 9 (*o,*o - 1)>; (19) 

r 2 = -^F^oK-; r 2j = ^^l(^(t)ln P<? (t ,to -*)>• (20) 

An unknown function p g (u) enters the expression f|T9l) through the terms In p(t) and 
In N. To get rid of this dependence, we use the averaging theorem. For the expressions 
for lnp(t) and lnA^ in (fT9l) we take all terms besides p q (u) out of the time integration. 
For each of these function however a different effective average time value should be 
used. The remaining integrals over p q {u) are equal to unity. We get: 

lniV~lnZi(c 3 ) -lnZi(c 4 ), (21) 
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\np(t)~]np q (c 1 ,c 1 -t) = -(52F m (c 1 )A m (c 1 -t)+]nZ(c 1 )). (22) 

in 

Let us make another approximation and change the order of the operations d/d(Aj) 

t 

and d/dt in the expressions (|T9l) - (!20l) . The value f r^dt = D(t) — D(t ) enters the 

to 

expression (fl8l) . where 



D(t) 



d 



((lnp(t) -lnJV) lnp 9 (t ,to-t)) 



= F i (to)[lnZ(c 3 ) - lnZ(c 4 ) - hiZ( Cl )] - 
Fi(ci) lnZ(to) + ^ {(A m A n ) - (A m )(A n )) 

m,n 

F m {Cl] 



(23) 



F n (to) 



(A i A n )-(A i ){A n ) 



_(A m Ai) — (A m )\Ai) 

(AkA m A n ) — (A m A n )(A k ) 



(AiA k ) - (Ai)(A k ) 



r 2j = — 



d_ 

Of 



rr, U » ' 



(y4fcAjA m ) — (AjA m ) (A k ) 



)-(A i )(A k ) 



+ lnZ(t c 



-E 



(^4jyl m ) — (A) (Am) 

(AjA m ) — (AjA m ) 



(24) 



The values of the correlators (Aj), (AiAk), (AjA k A m ) are averaged with p(t) and 
are t-dependent. In deriving (12"31 . f[2~4"|) we used the relations like 



dlnZ(c) 



E 



dlnZ(c) <9F n (c) 
9F n (c) 9(A) 



E 



(AA)-(A(A)' cj 



9(A) (A^ TO >-(A>(An)' 

One can proceed with the expression ( 1241) using the relations 



flfl _ ^ dFj d(A k ) _ v 8F k d(A k ) dF t 
dt ^ d(A k ) dt ^d(A) dt ' d(A k ) 

dlnZ _ _ sr^dFn 

dt ~ ^ dt [ mh 

m 

The solution to (fl8l) has the form 

p q (0)M 



dF k 



d 2 S 



d(A % ) d(A k )d(A) 



Pq(t ~ k 



1+^(0) J ^r 2 Mdt 



(25) 



So 
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M(t) = exp - t ) - 1 jf ndtj 



exp 



-§(«-*o)-^(I>(*)--D(«o)) 



(26) 



M(t ) = 1, where ridt is written in (123)) . 

Integrating by parts the second term in the denominator of fl25|) write it in the 
following form: 

p 9 (0)M(0. 



1 -L 



where 



/,**) -(/,**) 



M{t) 



to 



" K Fi dt 



(A Ac) - fWAfc) 



+ 



F;lnZ(t )-EE F ift 



(AyA m ) — (Aj)(A m ) 

{AiA m ) — (A) (An) 



(27) 
Af(*)dt 

(28) 



dD(t) 
dt 



E 



(A m A n ) — (A m )(A n ) 

(An A) — (Art) (A) 



eft 



ln( (A m A„) - (A m )(A n ))- 



p / \ (A m A n ) — (An) (At) 



-— ln((AnA) - (An) (A)) 
dt 

j t H(A m A n ) - (A m )(A n )) - lln((AA„) - (A)(A,)) 



-F m {ci)F n {t ) 



\AkA m A n ) — (A m A n )(A k ) 
{A l A k )-{A i ){A k ) 



d_ 

dt 



\ n ((AkA m A n ) ~~ (AnAi)(-<4 



lln((A^ fe >-(A)(^)) 



The expression for M(£) is given in (126|) . and p g (0) is determined from the conditions 
for the norm J Q °° p q (u)du = 1 and for the average lifetime (r) = J °° up q {u)du. 

If one either considers the stationary case or assumes a weak time dependence in 
the correlators in (j23j|-(j23l). D{t) ~ D(t ), r 2 ~ 0, and (125]l - (l27jl take on the form 



C 

P«(* - *o) = Pd(0) exp{— £(t - t )} 



(30) 



with Ci/Fi = p q (0) = l/(r). In ()3"U1) an exponential distribution for p q (u), used in 
[H [2J [3] is obtained. However generally the correlators in (I23]) - (I23| are time dependent. 
Applying the full form of (127)) for concrete systems it is possible to state the conditions 
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where the denominator in (I27p is considerably different from unity, hence the lifetime 
distribution essentially deviates from the exponential one f l30l) . 

The term L in the denominator of ( 1271) is small, since p q {0) ^ 1/ (r) <C 1. In the 
denominator of (f2"7|) it stands in the combination with 1, 1 — L « 1 . Therefore one 
can write p g = p g (0)M[l + L + L 2 + . . .] , and the distribution (1271) is close to the 
exponential distribution fl30|) used in [H [2l El [7J [8] . This results agrees with the results 
of [13], where the exponential distribution for the lifetime is shown to be a limiting 
distribution. But the expression for (r) is explicitely given in ( l30i) . and in ([211- (1291). It 
was already pointed above that the situations where L is comparable to 1 can arise as 
well. 

For the maximum entropy principle with the Shannon measure of the information 
entropy the exponential distribution used in [H [2], [3], is basic. Choosing another 
measures for the information entropy (e.g. Tsallis and Renyi measures [2J[25]) changes 
the function Jme ©, which yields another forms of the lifetime distributions. For the 
NSO method the functions Jme are in the form ([5]), and the information entropy is 
given by the Shannon measure, hence basic distribution is the exponential one. 

The expressions ( 1251) . ( 1271) for p q {u) depend on F i: the functions Ci = 
r k , k = 1,2, in (fl8|) - (fi9|) depend on the index i, the function p q (u,i) depends on i as 
well. One can get an ^-independent function p q (u) by symmetrizing the distribution, for 
example, using the operation 



Such formulation of maxent principle, as in [19J, gives the distribution for the 
lifetime, related to the exponential distribution which serves in this the base 

one. Distributions of other type can be obtained using some other form of maxent 
principle. 

3. Another approach to the method of maximum entropy 

Yet another approach to the determination of the type of function of distribution of 
lifetime is related to the method of maximal entropy inference ("maxent"), developed 
in [26] for the determination of the distribution of superstatistics. We note a formal 
similarity between the superstatistics method where the averaging is performed over 
the parameter (3, (for example, the inverse temperature) 



and the NSO method where the averaging is performed over the past life spans u = t—to, 



M 



Pq(u)= Y[p q (uri) 



i=l 
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which was already used in [27J. Therefore the maxent method of [2H] can be applied for 
determining the function p q (u). The analogy here is not merely formal. The principal 
assumption of [26j, that is the separation of the time scales is essential for the NSO 
method as well [H [2j [3] . In the approach of superstatistics [281 [23 EQ] the system is 
split into cells and local fluctuations of the value (3 are considered; the fluctuations of 
the value u — t — to affect the complete system. 

Closely related is also the research of Crooks [31] . He studies general non- 
equilibrium systems, without assuming that the system can be divided into different 
cells that are at local equilibrium. Crooks claims that instead of trying to obtain the 
probability distribution of the entire non-equilibrium system, one has to try to estimate 
the "metaprobability," the probability of the microstate distribution. Crooks also uses 
the maximum entropy principle but sets in ( !3~TT) A3 = 0. The main difference is that 
Crooks does not assume a local equilibrium in the cells, hence his approach, though 
being an interesting theoretical construction, does not give a straightforward physical 
interpretation to the fluctuating parameter (3. However such an approach can be applied 
to the fluctuations of the value u = t — t . In the approach of [26] one obtains a local 
fluctuating temperature that coincides with the thermodynamic temperature and which 
can in principle be measured. The work of Crooks is used by Naudts [32J to describe 
equilibrium systems. The author shows that some well-known results of the equilibrium 
statistical mechanics can be reformulated in a very general context with the use of the 
concepts introduced in [26], [29], [31] . 

In |26j following expression is obtained for the distribution function of the 
superstatistics f(/3; Aj) 

where Aj are Lagrange multipliers, V being an arbitrary constant (taking out a common 
factor out of the definition of Ai and A2 will turn out to be useful in the following). 
Using the well-known formula Sp = In Z(/3) + /3E(/3) with Z(Xi) being a normalization 
constant that is fixed by the condition (1)^ = 1. 

The same approach with said limitations used for the function p q (u), yields 

Z{t-u)-^' v E m Fm(t-u)(A m ) 
P q (u; Xi) = ex P( — A2 y A 3 </(u)), (31) 

where g(u) is an arbitrary function of u, which is determined by the physical peculiarities 
of the behaviour of the system in one or another period of its history. Expression fl3Tl) 
is obtained by the optimization of the entropy S(Xi) = J p q (u)\np q (u)du with the 
constraints for entropy 



p q (u)S(u)du = J p q {u) J p q (t — u,—u) In p q {t — u,—u)dzdu 
J P q (u)[-J2 F m(t-u)(A m )- In Z{t-u)\du 
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and parameters J p q (u) Yl m Fm{t—u)(A m )du. Similarly to [26], one can set the functions 
g(u) in a different fashion. For example, for g(u) = u, X 1 = A 2 = 0, A 3 = l/(r), 
where (T) is the average span of past life of a system (till the present time moment), 
we get the exponential distribution used in [TJ [2j [3]. Setting g{u) = ln-u with 
appropriate corresponding values of A one gets the power-like distribution for p q (u); 
setting g(u) = (Inu) 2 with corresponding A results in the log-normal distribution. Thus 
setting the function g{u) and A accordingly it is possible to obtain various distributions 
for the lifetime considered, for example, in |33j . 

It is possible to examine more difficult cases when the behaviour of the system 
changes at different stages of its evolution, when, for example the function g(u) = Inu 
yields the power-like function p q (u) at u < c, and g(u) gives an exponential shape of 
p q (u) at u > c. 

4. Conclusion 

For the determination of the lifetime distribution in the NSO method the method of 
maximum entropy principle as in [19] is used. The obtained distribution is close to 
exponential p qzu b(u) (T5]), but does not coincide with it. It is possible to find conditions 
at which this difference is essential. Using other variants of the maximum entropy 
principle, as in |26j, it is possible to obtain other distributions except exponential one, 
in particular, power-like and log-normal distributions, transitory ones between them, as 
well as distributions of other classes. 

In the interpretation of [2] it is the random value to in u = t — to that fluctuates. 
In [2] the limiting transition is performed for the parameter e, e — > in the exponential 
distribution p q (u) = eexp{— eu} after the thermodynamic limiting transition. In the 
interpretation of [TO] it corresponds to the average lifetime of a system tending to infinity: 
(T) = (i — to) = l/s —>■ oo. But the average intervals between successive random jumps 
grow infinitely, getting larger than the lifetime of a system. Therefore the source term 
in the Liouville equation turns to 0. If however the distribution p q (u) changes over 
the interval of the lifetime, the influence of the environment which caused this change, 
remains within the life span even if the lifetime tends to infinity [12]. 

There are numerous experimental confirmations for such change of the lifetime 
distribution p q (u) over the interval of the system lifetime. The examples thereto are 
the transition to chaos and the transition from laminar to turbulent flow which are 
accompanied by the change of the distribution of p q {u). In [MJ [35] the passage from 
Gaussian to non-Gaussian behaviour of the distribution of the first-passage time for 
some time moment is demonstrated. Besides the real systems possess finite sizes and 
finite lifetime. Therefore influence of surroundings on them is always present. 

Slow change of the function g{u) = u = t — t corresponds in the interpretation of 
[2] to the slow change of the random value to on a temporal scale. Accordingly slow 
change of the function g(u) = Inu corresponds to the slow change of ln(t — t ). Setting 
other functions g(u), for example, g{u) = (ln-u) 2 and so on is explained on the same 



Maximum entropy principle in NSO 



14 



footing. 

In the present work by means of two variants of the method of maximum entropy 
we obtained the expressions for the distribution of the lifetime value. It is noted that 
the choice of the form of the distribution function for the lifetime value can affect the 
non-equilibrium behavior of a system even after performing the thermodynamic limiting 
transition. 
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